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port  on  some  preliminary  experiments  in  the  use  of  this  code.  In 
particular,  we  have  computed  the  invariant  factors  of  all  (0,  1)  ma¬ 
trices  of  order  <  8,  with  constant  row  and  column  sums,  and  these 
data  are  presented  in  the  Appendix. 

An  examination  of  these  data  suggested  three  conjectures,  all  of 

♦  This  research  was  supported  in  part  by  the  Office  of  Naval  Research  under  Contract  No  Nonr  3775(00), 
NR  047040. 

NO  OTS 
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which  turned  out  to  be  true,  and  one  of  which  suggested  some  interesting 
questions  concerning  the  imbedding  of  a  non-singular  matrix  in  a  doubly 
stochastic  matrix  of  the  same  rank.  The  proofs  of  the  conjectures 
(Remark  1,  Remark  2  and  Theorem  2)  and  the  discussion  of  the  questions 
of  imbedding  (Theorem  1,  Theorem  3  and  Theorem  4)  form  the  main  part 
of  the  note.  We  hope  that  others  may  discern  additional  facts  from  the 
data  tabulated  in  the  Appendix. 

2.  Some  Simple  Remarks 

Let  Ot  be  the  class  consisting  of  all  m  x  n  (0,  1)  matrices 
with  prescribed  row  sums  and  column  sums. 

Remark  1.  The  ranks  of  matrices  in  LK  assume  all  integers 
between  the  minimum  rank  and  the  maximum  rank  of  matrices  in  CX  . 

Proof:  Let  A  be  a  matrix  in  CX.  Consider  the  2x2  sub¬ 
matrices  of  A  of  the  types 

Al“(o  l]  ind  A2=(j  j]  • 

An  interchange  is  a  transformation  of  the  elements  of  A  that  changes  a 
minor  of  type  A.  into  type  A  or  vice  versa  and  leaves  all  other 
elements  of  A  unaltered.  The  interchange  theorem  of  H.  J.  Ryser 
[  3  ]  states  that  if  A  and  B  belong  to  IX ,  then  A  is  transform¬ 

able  into  B  by  a  finite  sequence  of  interchanges.  We  note  that  if  B  is 
a  matrix  obtained  from  A  by  an  interchange,  then  B  =  A  +  C  where  C 


3. 


is  a  matrix  whose  entries  are  all  zero  except  for  a  minor  of  the  form 
,  and  therefore  rank  B<  rank  A+  1.  Similarly,  rank  A< 
rank  B  +  1.  Hence,  the  rank  is  altered  by  at  most  1  by  an  interchange, 
and  the  statement  to  be  proven  follows  at  once  by  the  interchange  theorem 
of  Ryser. 

From  now  on,  let  J  be  a  matrix  of  the  appropriate  size  whose 
entries  are  all  1,  and  let  A  =  J  -  A. 

Remark  2.  Let  A  beany  nxn  matrix  with  rational  integral 
entries.  Then  the  number  units  among  the  invariant  factors  of  A  and 
the  number  of  units  arm  ng  invariant  factors  of  A  differs  by  at  most 
1. 


Proof:  Let  A^  be  the  matrix  obtainable  from  A  by  subtracting 
the  first  column  from  cv. *ry  other  column,  and  let  *A^  be  the  matrix 

/V 

obtainable  from  A  corresp  u.in  .l’  Fxcept  possibly  for  the  first 

columns,  A^  and  are  the  same  j-  i  sign.  and  therefore  if 

e.,e_,,..,e  .with  e  lu,  |...|e  ,  a^e  the  invariai.,.  factors  of  A,  and 
12  n  1  1  2  11  n  1 

e’,  ,  e_, .  .  .  .“e  ,  with  I'e.  I. .  .  I'e  ,  are  the  invariant  fai  u  of  ,  then 
i  c.  n  1  c  n  1 


k 

n  e 
m=  1  1 


k+1 
n  "e 
m- 1 


and 


m 


k 

FI  e 
m=  1  r 


k-t  1 
n  e 
m=l  r 


k  -  1,2, 


Suppose  e.  is  the  first  non-unit  invariant  factor  of  A,  and  e.  is  the 
i  1  j 

first  non-unit  invariant  factor  of  A^.  We  may  as  well  assume  that 


i  <  j  and  j  >  2.  Taking  k  =  j  -  2,  we  have 


Cj-2 


j-2  i  j-1 

n  e  !  n  *?  =  +  1 

mi  m  - 

m=l  I  m= 1 


and  therefore  j  -  1  <  i  <  j,  which  completes  the  proof. 


3.  Imbedding  Questions 

Let  (X.  (  n,  k,  p)  denote  the  class  of  ail  n  x  n  matrices  with  real 
entries  whose  row  and  column  sums  are  all  k  ^  0,  and  whose  rank  is 
p  <  n.  We  shall  write  B  <  (X  (n,k,  p)  for  a  non-singular  matrix  B  of 
order  p  if  there  exists  a  matrix  A(  CX(n,  k,  p)  which  contains  B  as  a 
submatrix.  We  derive  as  Theorem  1  a  necessary  condition  for 
B<0((  n,k,  p).  As  a  corollary,  we  prove  a  relation  between  the  invariant 
factors  of  a  rational  integral  square  matrix  A  and  the  invariant  factors 
of  J  -  A.  In  Theorem  4,  we  show  that  this  necessary  condition  is  also 
sufficient  to  imbed  a  rational  integral  matrix  B  in  a  rational  integral 
matrix  A  e  CX(n,k,  p).  We  also  show  in  Theorem  3  a  set  of  necessary 
and  sufficient  conditions  for  imbedding  a  non-negative  matrix  B  in  a 
non-negative  matrix  A  e  CX(n,  l,p),  i.  e.  ,  a  doubly  stochastic  matrix  of 
rank  p. 

Theorem  1.  If  B  <  l/T(n,k,  p),  then  the  sum  of  the  elements  of 


/''B 

=  Vp  Ey 


Proof:  Assume 


5. 


where  A  e  t^|(n,  k,p).  Since  A  is  of  rank  p, 


(3.  1) 


E  =  DB  C. 


Let  be  the  column  vector  with  t  coordinates,  each  of  which  is  unity. 


Since  the  row  sums  of  A  are  k,  we  have 


(3.2) 


Bu  +  C  u  =ku, 
P  n’P  P 


(3.3) 


Du+Eu  =ku 
p  n-p  n-p 


Inserting  (3.  1)  and  (3.  2)  in  (3.  3),  we  obtain 


(3.4) 


D  B  u  =  u 

P  n'P 


Since  the  column  sums  of  A  are  k,  we  have 


(3.5) 


u'B+u  1  D: ku  1  . 
P  n-P  P 


Multiplying  both  sides  of  (3.  5)  on  the  right  by  B  u^,  and  substituting 


in  (3.  4),  we  obtain 


p+n-p=ku'B  u, 

P  P 


which  was  to  be  proved. 

Theorem  2.  Let  A  be  a  matrix  of  order  n,  with  row  and 

column  sums  k,  0  4  k  4  n>  whose  entries  are  rational  integers,  and  let 

A  =  J  -  A.  Let  "e . , . . .  /e  be  the  invariant  factors  of  A;  e, ,  .  .  .  ,  e 
In  In 

the  invariant  factors  of  A.  Then 


(3.6) 


n  e. 

cr<0  ‘  _  _k_ 

n  T  n-k 

T4o  1 


6. 


up  to  a  unit  -  1. 

Proof:  Let  A  be  of  rank  p,  B  a  non-singular  matrix  of  order 
p  contained  in  A,  and  B  =  J  -  B.  We  first  show  that  the  rank  of  B  is 
equal  to  the  rank  of  X  and  that  |B  |  4  0*  Write 

A  =  (d  Sj)  • 

Since  the  column  sums  of  A  are  all  k,  we  have  from  (3.  1) 

(3. 7)  u  1  C  +  u  1  D  B'1  C  =  k  u  \ 

p  n-p  n-p 

Substituting  from  (3.5),  we  have 

(3.  8)  u  ’  B_1  C  =  u 

P  n-P 

Let  X  be  a  matrix  with  p  rows  and  n-p  columns  such  that 

(3.9)  BX  =  C. 

From  (3.  8),  we  have 

(3.  10)  u  1  X  =  u 

P  n-P 

It  is  then  clear  from  (3.  9)  and  (3.  10)  that 
(3.  11)  BX  =  C. 

Further,  since  DX  =  E,  from  (3.  1),  it  follows  from  (3.  10)  that 
(3.12)  DX=E. 

In  other  words,  the  last  n-p  columns  of  X  are  linear  combinations  of 
the  first  p  columns.  Similarly,  we  can  see  that  the  last  n-p  rows  of 
X  are  linear  combinations  of  the  first  p  rows.  Consequently,  rank 
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=  rank  B  <  p  =  rank  B  =  rank  A.  But,  symmetrically,  rank  A  <  rank 
A,  which  implies  rank  A  =  rank  A  and  therefore  rank  H?  =  p,  from 
which  |1}j  0  follows.  Now,  we  have  from  Theorem  1  that 


(3.  13) 

We  shall  use  (3.  13)  to  prove  that 


.  _  - 1  n  ,  ~  - 1  n 

u'B  u  =  —  ,  u  1  B  u  =  — —  . 
p  p  k  p  p  n-k 


(3.  14) 


H  =  (.ijP-1  Jl_ 

B I 


Clearly,  (3.  14)  implies  (3.  6),  for  the  numerator  of  the  left-hand  side  of 
(3.  6)  is  the  g.  c.  d.  of  the  determinants  of  order  p  contained  in  A,  and 
the  denominator  is  the  g.  c.  d.  of  the  corresponding  determinants  in  3T. 
To  prove  (3.  14),  observe  first  that 


(3.  15) 


u  '  B  *  u  =  u'(B1u)  =  E  A., 
P  P  P  P  |B|  j  J 


where  A^  is  the  determinant  of  the  matrix  whose  kth  column,  k  j, 
is  B  ,  the  kth  columns  of  B,  and  whose  jth  column  is  u  .  Now, 

k  p 


(3.  16) 


=  V  'V  B2'B1 . V®!1’ 


which  can  be  verified  from  the  expansion  of  the  right-hand  side.  Further, 
if  we  apply  the  same  observations  to  *B,  we  have 


(3. 17) 


u  '  *§  *  u  -  ““57 


1 


|B|  J  4j  ' 


where  is  defined  in  the  obvious  way,  and 


8. 


(3.18) 


?  Aj  =  |up,  B2  -  Bj . Bp-  Bj|. 


But  Bk-Bj  *  -(B^-Bj),  and,  therefore,  from  (3.18)  and  (3.16),  we  have 

LA.  =  (-l)P-1  L  A.. 
i  3  j  } 


(3.19) 


Finally,  (3. 14)  follows  at  once  from  (3. 19),  (3.  17),  (3. 15)  and  (3. 13). 

Theorem  3.  Let  B  be  a  non- singular  matrix  of  order  p  whose 
entries  are  non- negative  real  numbers.  In  order  that  B  be  a  submatrix 
of  a  doubly  stochastic  matrix  A  of  order  n  and  rank  p,  where 
p  <  n  <  2p-l,  it  is  necessary  and  sufficient  that 


(3.20) 

(3.21) 

(3.22) 


u  'B“  *  u„  =  n, 


Sbij<l,  i  =  l,...,p;  L  b  jj  <  1 ,  j  =  !,...,  p, 


2  by  >  2  p  -  n. 


i.j 


Proof.  The  necessity  of  (3.20)  is  contained  in  Theorem  1.  The 
necessity  and  sufficiency  of  (3.21)  and  (3.22)  in  order  to  effect  the  imbed¬ 
ding  in  a  doubly  stochastic  matrix  without  considering  the  rank  was  point¬ 
ed  out  in  Dulmage  and  Mendelsohn  f  1  ] .  To  prove  the  sufficiency  of  (3.  20) 
-  (3.22)  with  the  rank  taken  into  consideration,  write 


A  =  (  D  E  )  • 


where  C,  D  and  E  are  to  be  determined.  Let  each  column  of  C  be 

-(I-B)u  ,  each  row  of  D  be  — -  u  '  (I—  B ),  and  each  entry  of  E  be 
P  p  .  P 


9. 


-  (  S  b..  +  n  -  2p).  Then  the  non -negativity  of  the  entries  of  A 

(n-p)2  i,j  1J 

follows  from  (3.  21)  and  (3.  22)  and  that  A  is  doubly  stochastic  is  easy  to 
verify.  To  show  that  the  rank  of  A  is  equal  to  the  rank  of  B,  we  must 
see  that 

— .  (u  '  B  u  +  n  -  2p)  =  -i-  u  '  (I-B)  B'1  -i-  (I-B)  u  , 
(n-p)2  p  p  "■>’  p  "-p  p 

which  follows  from  {3.  20). 

Theorem  4.  If  B  is  a  non- singular  matrix  of  order  p  with 
rational  integral  entries,  and  the  sum  of  the  coefficients  of  B  *  is 
^  ,  n  >  p,  then  B  is  contained  in  a  matrix  of  order  n  and  rank  p,  with 
rational  integral  entries,  whose  row  and  column  sums  are  all  k. 

Proof:  We  wish  to  find  a  matrix 


with  the  desired  properties.  Let  all  columns  of  C  other  than  the  last  be 
the  same  as  columns  of  B,  and  choose  the  last  column  in  such  a  way 
that  the  row  sums,  for  each  of  the  first  p  rows,  shall  be  k.  Next,  let 
all  rows  of  F  other  than  the  last  be  the  same  as  rows  of  (B,  C).  Choose 
the  last  row  of  F  in  such  a  way  that  all  column  sums  of  A  are  k.  All 
that  needs  to  be  checked  is  that  the  rank  of  A  is  p,  which  can  be  done 


as  in  the  previous  theorem. 
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APPENDIX 


How  to  read  the  table: 

(n,  k)  -►  (n,  n-k) 

81*  82’  • • •  •  Sn  hl*  h2 . hn 

means  that  if  A,  a  nxn  matrix  of  zeros  and  ones  whose  row  sums  and 

column  sums  are  all  k,  has  invariant  factors  8j»  82*  •  •  •  *  8n»  then  J-A 

has  invariant  factors  h, ,  h,,  . . . »  h  . 

x  L  n 

(n,  k)  -►  (n,  n-k) 


glg2' 


8 


n 


/hi,h  2 . 

\h'l ,  h  2  .  . . 


n 

,  h' 


. . n 

means  that  if  A  i  s  the  same  matrix  as  described  above,  then  J-A  has 


invariant  factors  either  h,,  h, . h  or  h',  ,  h',,  . . . ,  h'  . 

14  n  14  n 


(4.  1) 

— * 

V 

(4.  3) 

1 

1 

1 

1 

1 

1 

1 

1 

3 

(4.2) 

4 — 

1 

(4.2) 

1 

1 

0 

0 

1 

1 

1 

1 

0 

0 

1 

1 

1 

0 

I 

1 

1 

1 

0 

3. 


(8. 

1) 

f 

(8.7) 

1 

1 

1 

1 

1 

1 

1 

1  1 

1 

1 

1 

1 

1 

1 

1 

7 

(8. 

2) 

-u 

(8. 

6) 

1 

1 

1 

1 

0 

0 

0 

0 

1 

1 

1 

3 

0 

0 

0 

0 

1 

1 

1 

1 

1 

0 

0 

l 

0 

1 

1 

1 

1 

1 

3 

0 

0 

0 

1 

1 

1 

1 

1 

1 

0 

o  1 

1 

1 

1 

1 

1 

3 

0 

0 

1 

1 

1 

1 

1 

1 

1 

0  1 

1 

1 

1 

1 

1 

1 

3 

0 

1 

1 

1 

1 

1 

2 

2 

0 

1 

1 

1 

1 

1 

1 

1 

12 

0 

1 

1 

1 

1 

1 

1 

2 

2  1 

1 

1 

1 

1 

1 

1 

1 

12 

(8. 

3) 

(8. 

5) 

1 

1 

1 

1 

3 

0 

0 

» ! 

1 

1 

1 

1 

5 

0 

0 

0 

1 

1 

1 

1 

1 

3 

0 

0  1 

1 

1 

1 

1 

1 

1 

5 

0 

0 

1 

1 

1 

1 

1 

1 

3 

1 

0  1 

1 

1 

1 

1 

1 

1 

5 

0 

1 

1 

1 

1 

1 

1 

6 

0  1 

1 

1 

1 

1 

1 

1 

10 

0 

1 

1 

1 

1 

1 

1 

1 

1 

3 

| 

1 

1 

1 

1 

1 

1 

1 

5 

1 

1 

1 

1 

1 

1 

1 

6  1 

1 

1 

1 

1 

1 

1 

1 

10 

1 

1 

1 

1 

1 

1 

1 

15  1 

1 

1 

1 

1 

1 

1 

5 

5 

1 

1 

1 

1 

1 

1 

3 

3  ' 

1 

1 

1 

1 

1 

1 

1 

1 

15 

1 

1 

1 

1 

1 

3 

3 

3  J 

1 

1 

1 

1 

1 

1 

3 

15 

14. 
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